Spin Hall Effect of Double-Index Cylindrical Vector Beams in a Tight Focus

We investigate the spin angular momentum (SAM) of double-index cylindrical vector beams in tight focus. Such a set of beams is a generalization of the conventional cylindrical vector beams since the polarization order is different for the different transverse field components. Based on the Richards-Wolf theory, we obtain an expression for the SAM distribution and show that if the polarization orders are of different parity, then the spin Hall effect occurs in the tight focus, which is there are alternating areas with positive and negative spin angular momentum, despite linear polarization of the initial field. We also analyze the orbital angular momentum spectrum of all the components of the focused light field and determine the overwhelming angular harmonics. Neglecting the weak harmonics, we predict the SAM distribution and demonstrate the ability to generate the focal distribution where the areas with the positive and negative spin angular momentum reside on a ring and are alternating in pairs, or separated in different semicircles. Application areas of the obtained results are designing micromachines with optically driven elements.


Introduction
In micromachines, elements can be driven by light [1,2]. This requires designing optical tweezers appropriate for driven elements, depending on their shape, material, and motion trajectory. The work [3] discusses how to control the mechanical motions of various particles in optical tweezers under complicated actuation of optical forces and torques by tightly focused laser beams.
Typically, a light beam comes out from a laser with a Gaussian shape. Then, for certain applications, not only for optical trapping, but also for optical data transmission and laser welding, the beam should be converted to attain an on-demand shape. For this purpose, a huge branch of modern optics, laser beam shaping, is developed [4].
The beam shaping techniques are developed both within a resonator and outside it [5], by using refractional or diffractional optical elements. External beam shaping can be done to shape a beam that maintains this shape on propagation [6], or in some specific area, for instance, in the focal plane [7].
However, in optical trapping, there can also be a need not only to trap a particle at a certain point, but also to make her do some movements. For instance, the particle can be forced to travel along some trajectory, or rotate around its own center. Such a rotation occurs when the light possesses the spin angular momentum (SAM) [8], or nonlinear polarization. Thus, in addition to the task of shaping the beam intensity distribution, there can be a task of shaping the SAM distribution. In addition to optical trapping, the SAM can be used as information in optical data transmission [9].
In paraxial approximation, the intensity shaping can be done for a single transverse field component of a homogeneously polarized light. However, for shaping the SAM, both transverse components should be tailored with a controlled phase delay between them. The problem becomes more difficult under tight focusing conditions.
Recently, it has been observed that when a linearly polarized light beam is tightly focused, then, near the focus, areas occur with elliptic polarization [10]. Since the areas with negative and positive SAM are spatially separated, this phenomenon is a manifestation of the optical spin Hall effect. Later, the same effect was discovered for tightly focused high-order cylindrical vector beams [11]. In [11], the SAM is distributed mostly on a ring and consists of alternating areas with positive and negative values. The cylindrical vector beams have the Jones vector J = [cos mϕ, sin mϕ] with ϕ being the angular polar coordinate and m being the polarization order (for m = 1, radial polarization). A further generalization is a two-index polarization singularity with the Jones vector J = [cos mϕ, sin nϕ], where m = n [12], i.e., such a generalized vector field has different orders on the different axes. Recently, we studied such fields with V-points and for many values m and n we obtained the Poincare-Hopf index analytically [13].
In this paper, based on the Richards-Wolf theory [14], we study what happens with the SAM of a light field with the double-index polarization singularity in tight focus. We obtain an expression for the complex amplitude near the focus. Then, based on this expression, we derive the formula for the SAM and found that it can be nonzero only for the orders m and n of different parity. For analytical prediction of the SAM distribution, we decomposed the in-focus light field into the orbital angular momentum (OAM) spectrum and estimated the contribution of each angular harmonic. It turns out that if a light field being focused is not of ring shape and has a homogeneous or decaying Gaussian shape, then the OAM-spectrum consists mainly of mth and nth angular harmonics which exceed the other harmonics by an order of magnitude. It allows us to estimate the polar angles with zero SAM and thus to predict the SAM distribution. As an example, we obtained SAM distributions on a ring where the areas with positive and negative SAM occur in pairs.

A Light Field with a Double-Index Polarization Singularity near the Tight Focus
In [13], we investigated a generalization of cylindrical vector beams, when the polarization indices of the E x and E y field components were different. The amplitude of the electric vector of such a field is given by where E is the strength vector of the electric field, ϕ is the azimuthal angle in the source plane, (m, n) is the two-index polarization order, θ is the polar angle, describing the tilt of the light rays to the optical axis, A(θ) is the amplitude of the source field as a function of the axis tilt angle. Directions of the electric vectors are illustrated in Figure 1. In paraxial approximation, the intensity shaping can be done for a single transverse field component of a homogeneously polarized light. However, for shaping the SAM, both transverse components should be tailored with a controlled phase delay between them. The problem becomes more difficult under tight focusing conditions.
Recently, it has been observed that when a linearly polarized light beam is tightly focused, then, near the focus, areas occur with elliptic polarization [10]. Since the areas with negative and positive SAM are spatially separated, this phenomenon is a manifestation of the optical spin Hall effect. Later, the same effect was discovered for tightly focused high-order cylindrical vector beams [11]. In [11], the SAM is distributed mostly on a ring and consists of alternating areas with positive and negative values. The cylindrical vector beams have the Jones vector J = [cos mφ, sin mφ] with φ being the angular polar coordinate and m being the polarization order (for m = 1, radial polarization). A further generalization is a two-index polarization singularity with the Jones vector J = [cos mφ, sin nφ], where m ≠ n [12], i.e., such a generalized vector field has different orders on the different axes. Recently, we studied such fields with V-points and for many values m and n we obtained the Poincare-Hopf index analytically [13].
In this paper, based on the Richards-Wolf theory [14], we study what happens with the SAM of a light field with the double-index polarization singularity in tight focus. We obtain an expression for the complex amplitude near the focus. Then, based on this expression, we derive the formula for the SAM and found that it can be nonzero only for the orders m and n of different parity. For analytical prediction of the SAM distribution, we decomposed the in-focus light field into the orbital angular momentum (OAM) spectrum and estimated the contribution of each angular harmonic. It turns out that if a light field being focused is not of ring shape and has a homogeneous or decaying Gaussian shape, then the OAM-spectrum consists mainly of mth and nth angular harmonics which exceed the other harmonics by an order of magnitude. It allows us to estimate the polar angles with zero SAM and thus to predict the SAM distribution. As an example, we obtained SAM distributions on a ring where the areas with positive and negative SAM occur in pairs.

A Light Field with a Double-Index Polarization Singularity near the Tight Focus
In [13], we investigated a generalization of cylindrical vector beams, when the polarization indices of the Ex and Ey field components were different. The amplitude of the electric vector of such a field is given by where E is the strength vector of the electric field, φ is the azimuthal angle in the source plane, (m, n) is the two-index polarization order, θ is the polar angle, describing the tilt of the light rays to the optical axis, A(θ) is the amplitude of the source field as a function of the axis tilt angle. Directions of the electric vectors are illustrated in Figure 1. In our work [15], we have obtained expressions for the Cartesian components of a linearly polarized optical vortex, focused on an aplanatic system. If an optical vortex with a topological charge m is linearly polarized along the axis x, then, in the input plane, the electric field is given by whereas near the tight focus, the complex amplitude reads as where (ρ, ψ, z) are the cylindrical coordinates with the origin in the focus.
In Equation (3), the functions I ν,µ are defined as follows: where k = 2π/λ is the wavenumber of monochromatic light with the wavelength λ, f is the focal length of the focusing lens, α is the maximal tilt angle of rays to the optical axis, defining the numerical aperture of the aplanatic lens NA = sin(α), J ν is the νth-order Bessel function of the first kind. The same way, if such an input field is linearly polarized along the axis y, the field components near the tight focus are equal to The field with circularly symmetric amplitude distribution A(θ) and with polarization (1) can be decomposed into a superposition of four linearly polarized optical vortices: cos mϕ sin nϕ = 1 2 A(θ)e imϕ 1 0 + 1 2 A(θ)e −imϕ 1 0 Using this decomposition, we get the field components of the field with polarization (1) near the tight focus: The longitudinal component of the spin angular momentum is defined as [16]: For simplicity, we suppose that the functions A(θ) are real valued. Thus, all the integrals I ν,µ are also real valued in the focal plane (z = 0). Then, substituting the transverse field components from Equation (7), we get the following expression: This expression is cumbersome, but it reveals that the SAM is zero unless the polarization orders n and m are of different parity. This is in contrast with [13], where we obtain the polarization singularity index of light fields with polarization from Equation (1). This index is, vice versa, nonzero only for the orders n and m of the same parity.
When the parity is the same for both orders, then polarization is linear in the focus, since, according to Equation (7), both E x and E y are proportional to i m+1 (or i n+1 ), multiplied by some real-valued function. Near the center (r << λ), if n > m ≥ 2, the transverse components E x and E y are approximately proportional to the vector J = [cos (m−2)ϕ, -sin (m−2)ϕ]. If m > n ≥ 2, they are proportional to the vector J = [cos (n−2)ϕ, -sin (n−2)ϕ]. Thus, a saddle-type polarization singularity is generated in the center [17].
For n and m of different parity, simplifications yield This equation is hard to analyze without knowing which terms contribute the most. Thus, we need to decompose the near-focus field into the angular harmonics and study its OAM-spectrum.

Balance of Light Field Energy near the Tight Focus
At first, we study the OAM-spectrum of linearly polarized light after tight focusing. According to Equation (3), if the input field is polarized along the axis x, then the x-component of the electric vector near the focus consists of three angular harmonics, whose topological charges are m, m−2, m + 2. The y-component consists of only two angular harmonics with the topological charges are m−2 and m + 2, whereas the z-component also consists of two angular harmonics with the topological charges are m−1 and m + 1. Each harmonic is proportional to the function I ν,µ from Equation (4). Therefore, such a harmonic has the following energy W ν,µ (see Appendix A): This expression indicates that the energy of the angular harmonic is independent of the distance z from the focal plane and of the topological charge of the optical vortex that determines the index µ.
The integrals (11) can be evaluated analytically only in simple cases, but, nevertheless, the contribution of each angular harmonic can be estimated. For example, if the field being focused is a uniform field with a constant amplitude A(θ) ≡ 1, then According to Equation (3), the x-component E x consists of three angular harmonics with their amplitude being proportional to the functions I 0,m , I 2,m+2 /2, I 2,m-2 /2. The component E y is a superposition of harmonics described by the functions I 2,m+2 /2 and I 2,m-2 /2. Finally, the component E z is a superposition of harmonics described by the functions I 1,m+1 and I 1,m-1 . Therefore, the total energy of the focal field is The brackets in Equation (15) illustrate, respectively, the energies of angular harmonics in the E x , E y , and E z field components. Since W 2,m+2 = W 2,m-2 and W 1,m+1 = W 1,m-1 , we get Substitution of Equations (12)-(14) into Equation (16) yields After simplifications, we get This expression is exactly the size of a part of a sphere bound by the polar angle α. Thus, if a light field with unit amplitude is converging from a spherical surface with a numerical aperture sin(α), this field has exactly the energy given by Equation (18). This means that we found the balance when the energy of the input field is equal to the sum of energies of all angular harmonics of all three Cartesian components of the electric field in the focus.
In an extreme case, when the numerical aperture is close to the unit, i.e., α ≈ π/2, we get W 0,µ = (7/6)π f 2 , W 1,µ = (1/3)π f 2 , W 2,µ = (1/6)π f 2 . The whole energy is coinciding with the square of a hemisphere: W = W 0,m + W 2,m+2 + 2W 1,m+1 = 2πf 2 . Thus, the total energy W of the input field is distributed in the focal field in the proportions shown in Figure 2. One third W/3 goes into the longitudinal component E z , and 2W/3 goes into the transverse components E x (5W/8) and E y (W/24). The energy of the component E x is distributed into the mth-order angular harmonic (7W/12) and into the harmonics of the orders m − 2 and m + 2, each of the energy W/48. The energy of the component E y is distributed equally in the angular harmonics of the orders m − 2 and m + 2, each of the energy W/48. The energy of the component E z is distributed equally in the angular harmonics of the orders m − 1 and m + 1, each of the energy W/6.
Obviously, if the input field is linearly polarized along the axis y, the energy distribution is the same, but the main portion (7W/12, or 58%) goes into the mth-order of the angular harmonic of the y-component E y .
The above energy proportions can change if the field intensity is not homogeneous, i.e., if A(θ) = 1. However, if the amplitude function A(θ) decays from the center to the periphery, then the contribution of the mth-order angular harmonic becomes even greater. Indeed, for instance, if the aperture is bounded by an angle α, then the energy of the side angular harmonics of the orders m + 2 and m−2 relative to the energy of the central mth-order harmonic is W 2,m W 0,m = sin 6 (α/2) 1 − cos 6 (α/2) = sin 4 (α/2) 1 + cos 2 (α/2) + cos 4 (α/2) .
when α decreases from π/2 to 0, the numerator also decreases, while the denominator increases. Thus, this portion decays.
with the square of a hemisphere: W = W0,m + W2,m+2 + 2W1,m+1 = 2πf 2 . Thus, the total energy W of the input field is distributed in the focal field in the proportions shown in Figure 2.  Obviously, if the input field is linearly polarized along the axis y, the energy distribution is the same, but the main portion (7W/12, or 58%) goes into the mth-order of the angular harmonic of the y-component Ey.
The above energy proportions can change if the field intensity is not homogeneous, i.e., if A(θ) ≠ 1. However, if the amplitude function A(θ) decays from the center to the periphery, then the contribution of the mth-order angular harmonic becomes even greater. Indeed, for instance, if the aperture is bounded by an angle α, then the energy of the side angular harmonics of the orders m + 2 and m−2 relative to the energy of the central mth-order harmonic is when α decreases from π/2 to 0, the numerator also decreases, while the denominator increases. Thus, this portion decays. For instance, if α = π/2, then W 2,m /W 0,m = 1/7 ≈ 0.143, but even if sin α = 0.95, then W 2,m /W 0,m ≈ 0.057, i.e., almost all energy goes into the central mth-order harmonic. Now we apply the above introduced technique for obtaining the OAM-spectrum of a light field with the double-index polarization singularity (1). Actually, the input field consists of four optical vortices of the orders m, -m, n, and -n. In the tight focus, each of these vortices splits into the several harmonics with the above derived energies. We suppose here for simplicity that α = π/2 and that these harmonics do not coincide with each other. Thus, we have the energy distribution (OAM-spectra of the components E x , E y , and E z ), illustrated in Figure 3. Now we apply the above introduced technique for obtaining the OAM-spectrum of a light field with the double-index polarization singularity (1). Actually, the input field consists of four optical vortices of the orders m, -m, n, and -n. In the tight focus, each of these vortices splits into the several harmonics with the above derived energies. We suppose here for simplicity that α = π/2 and that these harmonics do not coincide with each other. Thus, we have the energy distribution (OAM-spectra of the components Ex, Ey, and Ez), illustrated in Figure 3. As seen in Figure 3, side angular harmonics of the orders m ± 2 and n ± 2 have relatively low energy (28 times lower than that of the orders m and n), and this energy becomes even lower when the aperture angle starts to decrease from α = π/2.

Spin Angular Momentum of Double-Index Polarization Vortices in a Tight Focus
Since it was found that almost all of the energy of the transverse field components goes into the mth-order and nth-order angular harmonics, we can suppose approximately that the SAM from Equation (10) reduces simply to Figure 3. Energy distribution of a tightly focused light field with double-index polarization singularity by the field components and by the angular harmonics. Numbers above each pillar indicate the fraction of the whole energy of the focused light field that goes to the given angular harmonic of the given field component (the sum of all numbers is 100%). Figure 3, side angular harmonics of the orders m ± 2 and n ± 2 have relatively low energy (28 times lower than that of the orders m and n), and this energy becomes even lower when the aperture angle starts to decrease from α = π/2.

Spin Angular Momentum of Double-Index Polarization Vortices in a Tight Focus
Since it was found that almost all of the energy of the transverse field components goes into the mth-order and nth-order angular harmonics, we can suppose approximately that the SAM from Equation (10) reduces simply to S z ≈ 2Im i n−m I 0,m I 0,n cos(mψ) sin(nψ).
Thus, it is seen that the SAM is equal to zero at the following polar angles: At certain conditions, these angles can coincide. For example, if n = 2m, we get S z ≈ 4Im{i m }I 0,m I 0,2m sin(mψ) cos 2 (mψ).
This indicates that there are 4m lines with zero SAM, starting from the origin and tilted with the polar angles with p = 0, 4m − 1. However, at odd p, the cosine in Equation (22) is zero, but it is squared. This means that the SAM does not change its sign at these angles. Instead, due to the square, there is a second-order edge dislocation. At even p, the edge dislocation has the first order, and the SAM changes its sign. In comparison with the first-order dislocations, the dislocations of the second order looks like wider dark areas between the maxima. Thus, the SAM distribution should look like a set of pairs of the spots with positive and negative SAM.
Another case occurs when m = 1. For the SAM to be nonzero, n should be even. Thus, the angles α = ±π/2 are again the lines of second-order edge dislocation, where the SAM does not change sign. At other angles, the SAM changes the sign.
If m and n are relatively large and close to each other, the roots of the sine and cosine do not coincide but are close to each other. Therefore, the SAM changes its sign at each such angle, but, due to the pairs of close zeros, the SAM between them is insignificant.

Simulation
Numerical simulation was done by the Richards-Wolf equations. At first, we studied the case when n = 2m. Figure 4 illustrates the intensity |E x | 2 + |E y | 2 + |E z | 2 and the longitudinal SAM density of a tightly focused light field with double-index polarization singularity of three different orders (m, n): (1, 2) (Figure 4a,d,g), (3,6) (Figure 4b,e,h), (7,14) ( Figure 4c,f,i) at the following parameters: wavelength λ = 532 nm, focal length of the lens f = 10 µm, numerical aperture sin α = 0.95, amplitude apodization function is homogeneous, i.e., A(θ) ≡ 1. The SAM density distribution was computed directly by Equation (8) and then it was compared with the one computed by Equation (10). The distributions were visually the same, relative error was computed as max|S z(10) -S z(8) |/max|S z(8) | (where S z (8) and S z (10) are, respectively, SAM densities computed by Equation (8) and by Equation (10)). Maximal relative error was at (m, n) = (7, 14) and equal = 1.9 × 10 −15 . The third row of Figure 4 illustrates the approximate SAM distributions obtained by Equation (20). Formally, the relative error from the distributions obtained directly by Equation (8) is large (14-32%), but it is seen in Figure 4 that it almost does not affect the shape of the distribution. This is because the error was caused by neglecting the side angular harmonics of the order m ± 2 and n ± 2, and since, according to Figure 3, there are 16 such harmonics, each of them separately is weak and cannot reshape the SAM distribution significantly. As seen in Figure 4, indeed, the SAM distribution consists of alternating pairs of spots with positive or negative SAM. This is different from the patterns we obtained earlier near the focus when the spots with positive and negative SAM were alternating in singles rather than in pairs [11]. Figure 5 illustrates the intensity and the longitudinal SAM density of tightly focused light field with double-index polarization singularity of two different orders (m, n): (6, 7) (Figure 5a,c,e) and (16,17) (Figure 5b,d,f) at the following parameters: wavelength λ = 532 nm, focal length of the lens f = 10 μm, numerical aperture sin α = 0.95, amplitude apodization function is homogeneous, i.e., A(θ) ≡ 1. For comparison, the SAM distributions were also computed directly by Equation (8) and approximately by Equation (20). The relative error is again large (14% for both fields), but the shape of the SAM distribution is almost undistorted. As seen in Figure 4, indeed, the SAM distribution consists of alternating pairs of spots with positive or negative SAM. This is different from the patterns we obtained earlier near the focus when the spots with positive and negative SAM were alternating in singles rather than in pairs [11]. Figure 5 illustrates the intensity and the longitudinal SAM density of tightly focused light field with double-index polarization singularity of two different orders (m, n): (6, 7) (Figure 5a,c,e) and (16,17) (Figure 5b,d,f) at the following parameters: wavelength λ = 532 nm, focal length of the lens f = 10 µm, numerical aperture sin α = 0.95, amplitude apodization function is homogeneous, i.e., A(θ) ≡ 1. For comparison, the SAM distributions were also computed directly by Equation (8) and approximately by Equation (20). The relative error is again large (14% for both fields), but the shape of the SAM distribution is almost undistorted. According to theoretical predictions, polar angles with zero SAM should occur by pairs of close angles. Figure 5 confirms it. It is seen that the positive SAM is mostly in the upper side while the negative SAM is mostly in the bottom side. Actually, the SAM is alternating, but between each spot with the positive or negative SAM, there is a weak spot of the opposite SAM, which is almost invisible in Figure 5.
The above theory predicts that the SAM is zero for the orders m and n of the same parity. Computation confirms this and polarization of the focal field is thus linear.  According to theoretical predictions, polar angles with zero SAM should occur by pairs of close angles. Figure 5 confirms it. It is seen that the positive SAM is mostly in the upper side while the negative SAM is mostly in the bottom side. Actually, the SAM is alternating, but between each spot with the positive or negative SAM, there is a weak spot of the opposite SAM, which is almost invisible in Figure 5.
The above theory predicts that the SAM is zero for the orders m and n of the same parity. Computation confirms this and polarization of the focal field is thus linear. It is seen in Figure 6 that Ey = 0 on the horizontal axis (φ = 0 and φ = π) and Ex = 0 on the vertical axis (φ = ±π/2), which is consistent with Equation (7) for the complex amplitudes of the light field. It is also seen that in both cases a saddle-type polarization [17] singularity is generated in the center.

Conclusions
Based on the Richards-Wolf theory, we have investigated here the spin angular momentum of double-index cylindrical vector beams in tight focus. Such a set of beams is a generalization of the conventional cylindrical vector beams since the polarization order is different for the different transverse field components. Thus, in the beam periphery, the number of areas with horizontal polarization is not equal to the number of areas with vertical polarization.
It turns out that if the polarization orders are of different parity, then the spin Hall effect occurs in the tight focus, which is there are alternating areas with positive and negative spin angular momentum, despite linear polarization of the initial light field.
On the contrary, if the polarization orders are of same parity, then polarization in the tight focus remains linear (but inhomogeneous).
For analytical description of the spin angular momentum distribution, we also analyzed the orbital angular momentum (OAM) spectrum of a linearly polarized mth-order vortex field in tight focus. It turns out that if the initial light field is not of ring shape and has a homogeneous or decaying Gaussian shape, then the energy of the angular harmonics with the orders m ± 2 in the transverse field components are at least 28 times lower than the energy of the mth-order angular harmonic.
This decomposition of the focused field into the OAM spectrum allowed us to predict the spin angular momentum distribution and, as an example, we demonstrated the ability to generate the focal distribution where the areas with the positive and negative spin angular momentum reside on a ring and are alternating in pairs or separated in different semicircles.
We limited our considerations by a homogeneous initial light field with constant amplitude and zero phase. Only polarization was supposed to be inhomogeneous. As a future work, light fields can be studied with the double-index polarization singularity, but with an inhomogeneous amplitude or/and phase distribution, like, for instance, circular Airy beams with vortices [18][19][20] or quadratic-power-exponent-phase vortex beams [21,22]. These beams were studied in homogeneous medium [18], uniaxial crystals [19,21], or in tight focusing conditions [20,22], but their tight focusing with double-index polarization singularities and the possibility of the optical spin Hall effect in the tight focus were not yet investigated. It is seen in Figure 6 that E y = 0 on the horizontal axis (ϕ = 0 and ϕ = π) and E x = 0 on the vertical axis (ϕ = ±π/2), which is consistent with Equation (7) for the complex amplitudes of the light field. It is also seen that in both cases a saddle-type polarization [17] singularity is generated in the center.

Conclusions
Based on the Richards-Wolf theory, we have investigated here the spin angular momentum of double-index cylindrical vector beams in tight focus. Such a set of beams is a generalization of the conventional cylindrical vector beams since the polarization order is different for the different transverse field components. Thus, in the beam periphery, the number of areas with horizontal polarization is not equal to the number of areas with vertical polarization.
It turns out that if the polarization orders are of different parity, then the spin Hall effect occurs in the tight focus, which is there are alternating areas with positive and negative spin angular momentum, despite linear polarization of the initial light field.
On the contrary, if the polarization orders are of same parity, then polarization in the tight focus remains linear (but inhomogeneous).
For analytical description of the spin angular momentum distribution, we also analyzed the orbital angular momentum (OAM) spectrum of a linearly polarized mth-order vortex field in tight focus. It turns out that if the initial light field is not of ring shape and has a homogeneous or decaying Gaussian shape, then the energy of the angular harmonics with the orders m ± 2 in the transverse field components are at least 28 times lower than the energy of the mth-order angular harmonic.
This decomposition of the focused field into the OAM spectrum allowed us to predict the spin angular momentum distribution and, as an example, we demonstrated the ability to generate the focal distribution where the areas with the positive and negative spin angular momentum reside on a ring and are alternating in pairs or separated in different semicircles.
We limited our considerations by a homogeneous initial light field with constant amplitude and zero phase. Only polarization was supposed to be inhomogeneous. As a future work, light fields can be studied with the double-index polarization singularity, but with an inhomogeneous amplitude or/and phase distribution, like, for instance, circular Airy beams with vortices [18][19][20] or quadratic-power-exponent-phase vortex beams [21,22]. These beams were studied in homogeneous medium [18], uniaxial crystals [19,21], or in tight focusing conditions [20,22], but their tight focusing with double-index polarization singularities and the possibility of the optical spin Hall effect in the tight focus were not yet investigated.
Application areas of the results obtained are designing micromachines for optical driving biological objects [23,24] or microtools in a lab-on-a-chip [25]. In contrast to the orbital angular momentum, which causes microscopic particles to rotate along a ring, the spin angular momentum causes articles to rotate around their centers of mass [8] and tailoring the SAM density distribution can allow simultaneous manipulation by an ensemble of particles. For this purpose, the particles should be non-metallic, but absorbing, since the intensity maximums attract dielectric particles. The particles should be of a size comparable with the areas of positive or negative SAM density, e.g., 0.2 µm in the distributions from Figure 4. Such sizes can have, for instance, polystyrene beads [26]. Another application is optical information transmission where the SAM density distribution can be used for encoding the data.

Acknowledgments:
We acknowledge the support of the Ministry of Science and Higher Education of the Russian Federation within a government project of the FSRC "Crystallography and Photonics" RAS (simulation).

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.